Quantum random walk : effect of quenching 
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We study the effect of quenching on a discrete quantum random walk by removing a detector 
placed at a position id abruptly at time tn from its path. The results show that this may lead to 
an enhancement of the occurrence probability at id provided the time of removal tn < t 1 ™ where 
t£ m scales as id . The ratio of the occurrence probabilities for a quenched walker (tn ^ 0) and free 
walker (tn — 0) shows that it scales as 1/tn at large values of tn independent of xd- On the other 
hand if tn is fixed this ratio varies as x 2 D for small xd- The results are compared to the classical 
case. We also calculate the correlations as functions of both time and position. 
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I. INTRODUCTION 



Quenching phenomena has been a much studied topic 
in both classical and quantum systems in the recent past. 
A quenching process can be broadly defined as one in 
which a certain quantity or condition of the system (e.g., 
temperature, magnetic field etc.) is changed in time 
from an initial value to a final value in a certain man- 
ner. Quenching can be studied in several possible ways. 
In certain cases, slow quenching is applied to obtain the 
classical ground states of the system as in the case of a 
spin glass system where there are many minima in the en- 
ergy landscape, separated by barriers which maybe over- 
come by quantum tunnelling [l|, On the other hand, 
in systems with a quantum critical point, nonequilibrium 
dynamics is studied by fast or slow quenching of relevant 
variables (usually a field or interaction which is made 
time dependent). Examples include transverse Ising and 
XY models 043 ■ I n ultracold atoms in an optical lattice 
fast quenching is applied by shifting the position of the 
trap potential and studying its response [f| 0] ■ In case 
of the transverse Ising or XY model, there is a deviation 
from the equilibrium state under quenching as the quan- 
tum critical point is crossed and the quantity of interest 
is the 'defect', which is the amount of departure from the 
actual equilibrium value. 



We consider here a discrete quantum random walk 
(QRW) la- tHl . completely different from a classical ran- 
dom walk [12h14| . on which the effect of quenching is 
studied. The QRW is made unitary by coupling the 
translation with chirality or rotation. The state of the 
walker is expressed in the |a;)|d) basis, where \x) is the 
position (in real space) eigenstate and \d) is the chirality 
eigen state (either left or right). 



II. QUENCHED QUANTUM WALK 
AND MEASUREMENTS 



SCHEME 



Usually, a quantum random walk is studied with two 
kinds of boundary conditions. In the infinite walk (IW), 
there is no boundary and in the semi-infinite walk (SIW) 
there is one absorbing boundary. Measurement-wise this 
signifies that there is no detector in the first case up to 
the time of observation, while there is a detector all the 
time in a particular position in the second. 

We consider the case when a detector initially placed 
at xd is removed suddenly at time This is termed a 
quenching phenomena as the presence of the detector and 
its subsequent removal may be compared to the presence 
of a time dependent transverse field in Ising model, XY 
model etc. We call this walk the quenched quantum walk 
(QQW). 

Classically, even when there is an absorbing boundary, 
a random walker will simply be either completely ab- 
sorbed there, or if not, will propagate like a free walker. 
In the quantum case on the other hand, the walk exists 
with non-zero probability at different locations and the 
absorption will thus take place with a probability. Ex- 
perimentally, this means that if the particle is detected, 
the evolution will be stopped. But if not, the walker will 
survive with non-zero and modified probabilities at the 
other locations [15| which is completely different from 
one which propagates without a detector [l6[ . 

In the QRW, the position of the particle, ip(x,t) is 



given as, 



i/j(x,t) = 



1pR.{x,t) 



(1) 



Here we have chosen the Hadamard coin [j| [l?} unitary 
operator H to perform the rotation which is coupled to 
the translation. H is given by 



H = J r 
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(2) 
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The occupation probability of site x at time t is given 
by |-0i(x,t)| 2 + \iPr(x, t) | 2 ; sum of these probabilities 
over all x is 1. The walk is initialized at the origin with 



2 




FIG. 1: (Color online) The probability of the walker getting 
detected at a point xd is shown against time. 



V>l(0,0) = a Q , iPr(0, 0) = 6 ; «o + ^o — 1- We have taken 
a symmetric walk; ao = 1/V%, bo = 

Let the detector be placed at some given site xd- In 
general, if the particle is at xd with probability a and 
the detector detects the particle with probability /3, then 
the total absorption probability at xd will be a/3. In our 
case we have chosen (3 = 1. Removing the detector at 
time tR is equivalent to having a step function behaviour 
of the probability of being detected at as a function 
of time as shown in Fig. Q] 

At each time step when the ensemble is measured the 
amplitudes describe only to the surviving copies. We 
may choose our scheme of measurement in two ways. Let 
the normalized occupation probability at x at time t be 
denoted by f(x,t). Thus f(x,t) denotes the fraction of 
the copies that survived the measurement up to time tn 
(not the fraction of the initial population) which reaches 
x at time t. If the particle survives absorption, then 
f(x,t) will be the conditional probability of finding it 
at x at time t in a 'single'observation. However, the 
average measure f(x, t) takes into account the absorption 
probability and is given by f(x,t)x survival probability; 

fi x t *) = 1 whereas J2x t) = l — d where d is the 
probability that it was absorbed earlier. This issue of 
two types of measurement was already addressed in [l6| . 

As xd and are the parameters of the system, 
we further modify our notation: f(x,t,XD,tii) and 
f(x, t,XD, to) are the normalised and average occupation 
probability of site x at time t respectively given xd and 
t R . 

Certain limiting cases can be immediately identified: 

(i) If tn = then the normalised and average occu- 
pation probabilities become same and are identical to 
the usual occupation probability of a quantum random 
walker when there is no detector at all; the IW case. 

(ii) If tn = oo, it is a case of SIW. 

(hi) If xd >> 0, the walk will be IW-like in finite 
times. 

The study for quench must be for a time t > tu > xd, 
otherwise the removal of the detector does not affect the 
probabilities. [Note that, here the numerical value of xd 
is considered while comparing with time t.] 
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FIG. 2: (Color online) Snapshots of the probabilities of occu- 
pation / for a symmetric walker under different conditions are 
shown for t = 100 (left panel) and 500 (right panel) against 
position x. (a) and (b) are results for IW. In (c) and (d) 
the results are shown for a SIW when a detector is placed 
at id = 10 for all times. In (e) and (f) QQW results are 
shown where the detector is initially placed at xd = 10 but 
is removed at time ta — 50. 



III. RESULTS 

We first present some snapshots of the probabilities of 
occupation f(x) of different sites at different times for 
the IW, SIW and QQW in Fig. [5] to get a comparative 
picture of the dynamics in the three cases. For a quantum 
walker, the displacement in time t is proportional to t and 
the maximum of the occupation probability occurs at a 
value of x ~ t/ y/2. These features are clearly shown for 
the IW. For the SIW and QQW, the detector placed at 
Xd will make the probabilities different; up to t = tn, 
QQW and SIW are equivalent. However, for t > tji, as 
shown in the figure, the probability "spills out" beyond 
xd for the QQW, although far away from xd, there is 
little difference. In fact away from the boundary, we find 
that maximum probability is again at a value of \x\ ~ 
t/V2. 

We next present the data for a fixed value of x = 
xd at different times t for both f(x,t,XD,tn) and 
f(x,t,XD,tn), using the shorthand notation / and / for 
these two quantities respectively. To study the effect of 
quenching, the ratios of the probabilities at different po- 
sitions and times for the IW and the QQW may be cal- 
culated for given values of xd and t^. This is in tune 
with the measure of defects in studies of quenching in 
the quantum spin models. The ratios f/fo and ///o are 
computed where /o is the occupation probability for the 
IW. f/fo shows that it can attain values much larger 
than 1 for a short time t > tu before saturating at larger 
times to values larger or smaller than 1 depending on tR. 
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FIG. 3: (Color online) Top: Ratio of f(xD,t,XD,tii) and 
fo for x D = 10 and t R = 20,30,50,100. Bottom: Ratio of 
f(xD,t,XD,tn) and fo for xo = 10 for the same values of xd 
and tR. 
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FIG. 4: (Color online) Plots of (f/fo)sat against ta for xo = 
10, 15, 20, 25 respectively. For large tjj the log-log plot shows 
a slope of —1. 
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These are shown in Fig. [3J 

The ratio of f/fo are also clearly different from 1 be- 
yond t — tR. However, it shows a much more regular 
behaviour compared to f/fo although it takes a longer 
time to saturate (Fig. [3]). As the detector is placed at 
xd, up to t = xd, the walker is not affected by the pres- 
ence of the detector and the ratios are equal to 1 for all 
tu. Beyond xd < t < tR, the initial population of the 
ensemble decreases due to detection at site xd- As time 
goes on, larger fraction of the initial population get lost 
and the ratio thus decreases gradually. Then as the de- 
tector is removed, the probability grows at xd, such that 
the ratio starts increasing before reaching the saturation 
value. In fact the results are of interest mainly at t > tR, 
after the detector is removed, when the QQW is different 
from the SIW and the ratio starts growing. Interestingly, 
contrary to naive expectation, the ratio may saturate at 
values larger than unity for small tR. The explanation for 
(/ / fo)sat having value > 1 is given later in the paper after 
we present the results for x ^ xd ■ The saturation values 
of the ratio, (f/fo) sat plotted against tR shows an ini- 
tial non-monotonic behaviour, but for large tR, (f / fo)sat 
clearly scales as 1/tR for all values of xd (Fig. 0J. 

From Fig. 21 we note that there is a value of tR beyond 
which (f / fo)sat is never > 1, i.e., beyond this particu- 
lar tR = t R m , (f / fo)sat monotonically decreases below 1. 



We find that t 1 ™ 



varies as x D (Fig. [SJ. Moreover, for 



fixed tR, the saturation value of the ratios show a vari- 
ation with xd- for large tR, (f / fo)sat varies as xd 2 , for 
smaller values of tR, this variation is valid over a small 
range of xd- Combining the above two results, we con- 



FIG. 5: (Color online) Plot of t| m against x D . The log- log 
plot shows a slope of 2. In the inset (f / fa)sat against xd is 
plotted for different ta. 



elude, for large tR, 



(f/fo) sat = kx D 2 /tl 



(3) 



where A: is a constant with dimension of inverse length. 
That (f/fo)sat should decrease with ta is expected as 
the probability of detection increases with larger tR. On 
the other hand as xd is increased, the walk remains un- 
affected for longer times which implies that (f/fo) sat 
should increase with xd (xd — > oo makes / = fo). How- 
ever, the exact scaling form eq. ([3]) is not obvious. 

Before discussing other issues it is worthwhile to com- 
pare the results of the quantum case to the classi- 
cal quenched random walk (CQW) under identical fast 
quenching. For a classical random walker, the occupation 
probability f c (x,t) is simply given by f c o (x, t)x survival 
probability, where f c o(x,t) is the occupation probability 
for the classical walker in absence of any boundary. The 
survival probability is given by l — J Q tR F c (xD,t)dt, where 
F c (xD,t) is the first passage probability at xd at time t. 
So for the classical case, 



A 

fcO 



F c (x D ,t)dt, 



(4) 



and is always less than 1. Moreover, it is independent of 
x and t. fdfco for tR — > oo scales as A= [HI]. For the 
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FIG. 6: (Color online) Ratio of the occupation probabil- 
ities f(xD + r)/fo(xD + r) against r for xd = 10 and 
tn — 20, 50, 100 at a particular time t = 100. Inset shows 
g/go against t for tR = 50 for different r. 

quantum case in contrast (/ / fo) S at has a different scaling 
behaviour with xd/V^R as already shown. 

In fact for a CQW, the ratio is simply identical to the 
persistence probability at time tR, i.e., the probability 
that the walker has not visited the site xd till time tR. 
For the quantum walker, the corresponding persistence 
probability is proportional to [lg]. Hence the satu- 
ration ratio is not linearly dependent on the persistence 
probability in case of the QQW in contrast to the CQW. 

Let us now discuss the results for sites x ^ xd) for a 
site separated by a distance r from xd the occupation 
probability can be written as f[xn + r,t, xd^r). The 
ratio of the occupation probabilities J(xd + r)/ fo{xD + 
r) (in short hand notation) shows remarkably different 
behaviour for r > and r < at a finite value of t (Fig. 
[6]). For r > 0, the ratio goes to zero at a finite value 
of r, the decay is smooth for large values of tR while for 
small tR, the decay is accompanied by small oscillations. 
For r < 0, the ratio shows an irregular behaviour with 
r, several peaks occur with the peak values much greater 
than 1. However, looking carefully at Fig. [21 it is evident 
that for sufficiently large \r\, when r < 0, the QQW and 
SIW behave in the same way, and the ratio is not much 
affected by the removal of the boundary. Hence one can 
say that the memory effects are strong here. 

As function of time, the product of {^ D+ , r \ and {^ n \ 
are plotted in the inset of Fig. [5] for xd = 10 and tn = 
50. The product approaches unity for all r, although 
for r > it reaches unity from below and for r < 
from above, this product can be identified as the ratio of 
two correlation functions g(xD+r,t,XD, tR.) and goixo + 
r, t, x u , tij) where 

g(x D + r, t, x Dl t R ) = f(x D + r, t, x D ,t R )f(x D , t, x D ,t R ) 

(5) 

and go(xn + r, t, XD,tR) is obtained by replacing / by /q 
in eq.©. Thus we note that the ratio of the correlation 
functions approach unity irrespective of the value of r for 
finite values of tR. 

The ratio of the two correlation functions may also be 
estimated for the classical case. Using the result of eq.(j4|), 



jj^ is simply the square of the persistence probability 
and is less than 1 always irrespective of the value of r 
and t. This is clearly different from the quantum case 
where the ratio goes to 1 for all r at large t. Moreover, 
in the quantum case there is a time dependence which is 
strongly dependent on the sign of r at small t. 



IV. SUMMARY AND DISCUSSIONS 

In summary, we have studied quenching in a quantum 
system in a completely different sense compared to earlier 
works. The observation that the occurrence probability 
of a QRW may actually be enhanced by quenching is one 
of the main results of this study. This is a purely quan- 
tum mechanical effect. Having the detector up to time 
tR means the quenched walker cannot go beyond x d , and 
the undetected walker will move away from x d ■ However, 
at a later time t > tR, when the walker is free once again, 
it can move towards xd and go beyond. From Fig. [2j it 
can also be seen that that most of the contributions to 
Xd and beyond come from the density of walkers closer 
to it, as the occurrence probability far way from xd are 
not much affected with the removal of the detector. At 
later times after the removal of the detector, the occu- 
pation probability profile approaches the IW picture as 
the local hill like structures smoothen out. This happens 
closest to at earlier times and slowly the further parts 
are affected. In comparison, in the infinite walk case, the 
walker has moved reasonably away from xd at tR, such 
that the ratio can be greater than unity close to xd- But 
when tR is greater than tR hm , the ratio can no longer 
exceed unity. At even later times, the ratio saturates as 
the "memory" of the detector gets erased in time. On 
the other hand we find that memory effects are strong 
for x < < xd where the removal of the detector is more 
or less irrelevant. Hence the effect of quenching is rather 
local. 

Other important results are the scaling behaviour of 
the quantities like tR hm with xd and (f/fo)sat at xd as 
a function of xd or tR. Although for QRW, the displace- 
ment varies linearly with time, we find that the timescale 
t R Um var j es w ith in a quadratic manner. The scaling 
behaviour of (f/fo)sat is also drastically different from 
the classical case. 

The present work can be extended in many ways like 
making the probability of detection dependent on time in 
a different way and slow quenching can be studied when 
it decays algebraically. One can also, instead of having a 
detector from time t = to t = tR, manipulate both the 
times at which the detector is placed and subsequently 
removed. Quantum random walks of correlated particles 
have been shown to encode information [l8j and quantum 
walk is capable of universal quantum computation [l9| . 
Quenching of quantum random walks may lead to some 
new features in such contexts, although it is too early to 
predict exactly how. 
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